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During eternal inflation, a landscape of vacua can be populated by the nucleation of bubbles. 

These bubbles inevitably collide, and collisions sometimes displace the field into a new minimum in 

a process known as a classical transition. In this paper, we examine some new features of classical 

transitions that arise when gravitational effects are included. Using the junction condition formalism, 

we study the conditions for energy conservation in detail, and solve explicitly for the types of allowed 

classical transition geometries. We show that the repulsive nature of domain walls, and the de Sitter 

expansion associated with a positive energy minimum, can allow for classical transitions to vacua 

of higher energy than that of the colliding bubbles. Transitions can be made out of negative or 

t_^ zero energy (terminal) vacua to a de Sitter phase, re-starting eternal infiation, and populating new 

^~~^ vacua. However, the classical transition cannot produce vacua with energy higher than the original 

^^ parent vacuum, which agrees with previous results on the construction of pockets of false vacuum. 

Cn We briefly comment on the possible implications of these results for various measure proposals in 

^>^ eternal inflation. 
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I. INTRODUCTION 

The ubiquity of vacua in theories beyond the standard model of particle physics has highlighted the importance 
of understanding how our particular low-energy vacuum was accessed in the early universe. The need for such an 
understanding is particularly urgent in theories with extra dimensions, such as string theory, where vast collections 
of low-energy vacua (features in a potential "landscape" ) can result upon compactification. A central question is to 
determine how (or if) the various states in the landscape can be populated. 

Typically, the vacua in a potential landscape will not be completely stable. Their decay occurs via the nucleation 
of bubbles containing some new phase, which subsequently expand into the parent vacuum. The process of bubble 
nucleation can in many cases be described by the Coleman de Luccia (CDL) instanton T-'3' . Starting in some positive 
energy minimum, the resulting inflationary expansion can end locally due to bubble nucleation events. However, when 
the rate of bubble formation is roughly less than one bubble per Hubble volume per Hubble time, the phase transition 
cannot complete, and inflation becomes eternal (for some reviews of eternal inflation, see [HIS]). Bubble nucleation 
during eternal inflation can populate many of the vacua in a landscape. 

Including gravitational effects in the description of domain walls and vacuum bubbles leads to some dramatically new 
phenomena. One consequence is that domain walls become repulsive [Si [7]. Another important, purely gravitational, 
effect is the existence of surviving bubbles containing false vacuum (a higher energy phase than the background). In 
the absence of gravity, the pressure gradient pointing into a false vacuum bubble will cause it to collapse to zero size. 
Including gravity, if the interior of the false vacuum bubble has positive energy, and is larger than the associated de 
Sitter horizon size, it cannot (by causality) collapse to zero size. 

Although the phase transition does not percolate, there will still be collisions between bubbles. In fact, since each 
bubble grows to reach infinite size, all bubbles will undergo an infinite number of collisions. Instead of dissipating 
the collision energy, the field configuration can remain compact and naturally form a domain wall enclosing a region 
of some other vacuum. Locally, the collision has kicked the field into a new basin of attraction. This represents 
a secondary mechanism for populating vacua during eternal inflation, and is referred to as a classical transition [S]. 
Classical transitions have been studied in simulations without gravity [8. , and their field dynamics in fiat space further 
elaborated on in [5]. 

These studies of classical transitions in flat space suggest that, in the thin- wall limit, the spacetime can be modeled 
as vacuum regions connected across domain walls. Under this approximation, we can apply the Israel Junction 
condition formalism [TU] to determine the global structure of the spacetime, including gravitational effects. The 
junction condition formalism has been an important tool in determining the outcome of bubble collisions and the 
large-scale structure of the resulting spacetimes. Some representative examples in the literature include |11H15| . 

In this paper, we explore the large-scale structure of classical transition geometries, determining the general criteria 
that must be met for a classical transition to occur, and outlining solutions where the effects of gravity are important. 
There are many parallels with the types of new solutions that arise for vacuum bubbles. One particularly interesting 
result is that classical transitions can give rise to a lasting region with a vacuum of higher energy than the interior of 
the colliding bubbles. This allows for the creation of inflating vacua from the collision of bubbles containing so-called 
terminal vacua, thereby providing a means of connecting otherwise isolated sectors of a landscape. We formulate 
constraints on the types of classical transitions that are allowed by energy conservation, and assess the types of vacua 
it is possible to populate using this mechanism. 

The plan of the paper is as follows. In Sec. [H] we briefly describe the field dynamics of classical transitions. We 
then describe how the Israel Junction condition formalism can be used to construct classical transition spacetimes 
in Sec. |III[ In Sec. |IV| we specialize to the collision of bubbles containing a positive or zero energy phase, classifying 
the possible classical transition geometries, and deriving a number of constraints from energy conservation. Some 
of these results are generalized and expanded upon in Appendix [A] and |B] Colliding bubbles containing a negative 
energy phase allows for a wider variety of classical transition geometries, which we study separately in Sec. IV] Some 
technical results are collected in Appendix [C] Finally, we briefly discuss the implications of these results and conclude 
in Sec. ED 

II. CLASSICAL TRANSITIONS: FIELD DYNAMICS 

Analytic studies of bubble collisions P^HTT] depend heavily on the assumptions made about how colliding domain 
walls interact. Some motivated guesses for the possible outcome of a collision between two bubbles include the 
creation of shells of radiation or dust, and in the case of non-identical bubbles, a domain wall that separates the 
bubble interiors. Numerical simulations of bubble collisions in fiat space [11] [18] have confirmed these types of 
behavior, but also revealed a number of new possibilities. Simulations by Easther et al [S] uncovered the existence 
of a phenomenon known as the "classical transition," whose study will be the main focus of this paper. Progress in 




FIG. 1. Left: A classical transition. Vacuum bubbles Bi and B2 collide and create a new vacuum region C. 

Right: A scalar field potential that can give rise to a classical transition with Bi — B2 — B. Each of the three minima is 

separated by a distance S(l)/2 in field space. 
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FIG. 2. Field configurations before(top) and after(bottom) a classical transition. The solitonic profiles superpose, and because 
of the equal spacing of the vacua in the potential, create a region in the C phase which subsequently expands. 



understanding the scalar field dynamics, presented in a recent paper [5], will allow us to set up our analytic study. 
We quickly review the main results of this study here. 

As shown in Fig.[l} a classical transition refers to a collision that creates two new domain walls enclosing a region of 
some new vacuum. This was already observed in early simulations by Hawking, Stewart, and Moss jlT]. who studied 
a potential with two minima, and found that a collision can cause the field to jump back up to the false vacuum in 
a localized region. With a more complicated potential, containing more than two minima, the there are many more 
possibilities. 

The prototype of the classical transition we study is realized by the potential in Fig. Ill which has three equally 
spaced vacua A, B, and C. Fig. [2] shows a numerical simulation of the prototype classical transition. Surprisingly, 
when the domain walls separating vacuum A from vacuum B collide, the field profiles can stay compact. When the 
velocity of the walls is large enough, the potential for the field becomes irrelevant as compared to the spatial gradients, 
and the solitonic field configurations simply superpose. Because the vacua are equally spaced, the domain walls now 
enclose a region of vacuum C, which subsequently expands due to its lower energy density as compared to vacuum BI 
It is also clear from the fact that the walls simply superpose that collisions can only cause excursions in field space of 
a fixed distance. 

In a classical transition, the resultant field profile may not exactly be the soliton between (j^g and (pc^ because the 
barriers are different. In particular, the height difference implies that the solitons have different mass (wall tension). 
In high resolution simulations [H], it was observed that the profiles settle into the true solitonic configuration, in 
the process radiating very little energy. Instead, the velocity is adjusted to accommodate the mass change. The 
superposition condition, which comes from microphysics, is always more restrictive than energy conservation, although 
both give a lower bound on the incoming boost necessary to produce a classical transition. 

In general, the classical transition in Fig. Ill can be realized on a potential for multiple variables (/) if the following 
conditions are true. 



• The incoming boost is large enough so that the soUtonic profiles freely pass through each other. 

• (f)Q = cf)g^ + 0^2 — 0^ is a vacuum. 

• The solitonic profile between A and Bi can settle into the new profile separating B2 and C. Similarly A and 
B2 to Bi and C. 

The last two conditions may look finely tuned, but they can be easily arranged. [J In the remainder of the paper, we 
will have little to add to the discussion of how ubiquitous classical transitions are in a realistic potential landscape. 

III. CLASSICAL TRANSITIONS WITH GRAVITY 

Neglecting gravity, the vacuum energy on either side of a domain wall fully determines its dynamics. In a spacetime 
undergoing a classical transition, as depicted in Fig. IT] the region of new vacuum (vacuum C) will expand indefinitely 
when its energy is lower than that of the colliding bubbles (vacuum B). If the region of new vacuum has higher 
energy, it is always doomed to collapse to zero size. Therefore, neglecting gravity, classical transitions can only seed 
lasting regions of new vacuum whose energy is lower than the energy inside of the colliding bubbles. 

How might the inclusion of gravity affect this picture? A full treatment of gravitational effects in a classical 
transition would require simulations employing numerical relativity. However, given a scalar potential with the 
appropriate vacuum structure, the field dynamics outlined in the previous section (more details will be presented 
in [9]) strongly suggest that it is reasonable to think of the classical transition spacetime as consisting of regions of 
vacuum separated by domain walls. We can therefore make an important simplification, and reduce the full field 
dynamics to a study of the motion of domain walls in a background spacetime. This simplified picture allows us to 
isolate several distinct phenomena: the gravitationally repulsive nature of domain walls [HI 13 1 the effects of a non-zero 
cosmological constant on the motion of domain walls, and the gravitational backreaction of the energy released in the 
collision. 

As we will see, the presence of gravity alters the domain wall motion in a nontrivial way, modifying the energy 
conservation condition q when the walls cross each other. However, the explicit form of energy conservation and 
the equations of motion is known and after working though this, the nontriviality in some cases conspires to allow 
for classical transitions to expanding regions with vacuum energy higher than that of the colliding bubbles. Such 
transitions are particularly dramatic when the colliding bubbles have negative vacuum energy, since the region of new 
vacuum can avoid the big-crunch singularity inside the bubbles, bringing a "terminal" vacuum back to life. 

In the remainder of the paper, we will examine these features using the Israel Junction condition formalism. After 
describing the general setup and the relevant equations of motion and energy conservation in this section, we move 
on to study colliding bubbles with positive and negative vacuum energy in Sec. |IV| and [V| respectively. 

A. General setup 

The collision between two vacuum bubbles preserves an S0(2,l) symmetry [TT1[T2]. Using a hyperbolic version of 
Birkhoff's theorem, the most general vacuum solution to Einstein's equations can be written as 

ds'^ = -ai{zy^dz^ + ai{z)dx^ + z^dHl , (1) 

where dH2 = dy^+smh X'^'^^ is the metric on a unit hyperboloid, and the coordinates generally range from < z < cx), 
—(30 < X < 00, —00 < X < 00, and < </) < 27r. Allowing for vacuum energy, the metric functions ai{z) are 

ai{z) = l-'^-^+HJz\ (2) 

z 

where the label / keeps track of different regions. In our conventions, _ff| > for positive vacuum energy, and Hj < 
for negative vacuum energy. 

The setup for the collision problem is shown in Fig. [31 which depicts a slice in the x, z plane. In the body of the text, 
we consider collisions between identical bubbles only, and briefly comment on the collision between different bubbles 



^ In a toy model with fluxes and extra dimensions llQI , these conditions are guaranteed. 

^ Energy conservation alone is not enough to guarantee a classical transition, since the conditions on the field dynamics can be much 
more stringent. In this paper, we will assume that a classical transition can proceed as long as the energy conservation is satisfied. 
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FIG. 3. The potential (left) and collision spacetime (right) that we consider for a classical transition. Two bubbles of the B 
phase (of initial radius 7?o) in a background of the A phase collide at z = Zc, yielding a classical transition to vacuum C. The 
C phase can have either higher or lower energy than the B phase. There might be a small amount of additional debris released 
by the collision, which we account for as null shells, enclosing region B. Each region, described by the metric Eq. ([l|, is sewn 
together using the junction condition formalism. Our goal is to find the fate of region C. 

in Appendix |B] We assume that the underlying scalar potential has three vacua, A, B, and C and can support a 
classical transition in the absence of gravity from vacuum B to vacuum C. We take H^ > Hg so that true vacuum 
bubbles of the B phase, with some initial radius Rq ^ H^^, can nucleate from the A phase. When two bubbles 
of the B phase collide at z — Zc (in these coordinates, the collision occurs at a position specified by z), a classical 
transition can occur to the C phase, which we allow to be of arbitrary energy H^. A domain wall of normalized 
tension kBC (this quantity is related to the tension through the relation ksc — ^T^Gasc) separates the B and C 
phases. In addition to the classical transition, the collision can release energy in other forms, which we model as null 
shells of radiation. We refer to the portion of the B phase to the future of a radiation shell as the B phase. While we 
generally include their effects, we neglect these energy sinks in many of the specific examples to follow. 

The metric in each of the four qualitatively different regions is of the form Eq. ([I]) , with the metric coefficient in 
each region given by 
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There are two free parameters in the metric, Mc and Mg, that are not determined by the underlying potential. 

The junction condition formalism gives a prescription for sewing each of these spacetime regions together in a 
consistent way. First, we require that the metric be continuous across all junctions. This implies that z must be 
continuous throughout the collision spacetime. The derivatives of the metric can be discontinuous, and will be related 
to the energy density of local sources such as domain walls. In the collision spacetime, there are a number of separate 
interfaces to consider: the bubble walls between the A and B phases, the post-collision domain walls separating the 
B and C phases, and the radiation shells separating the B and B phases. The dynamics of each interface can be 
determined from Einstein's equations. Requiring that each of these interfaces are stitched together consistently is 
equivalent to imposing energy and momentum conservation at the location of the collision. Summarizing the necessary 
steps in solving for the collision spacetime, we must: 

1. Solve for the trajectories of the in-coming bubble walls and the location of the collision. We work in a frame 
where the colliding bubbles are nucleated at the same time, and so the kinematics of the collision are fully 
determined by their initial radii and separation (see jl4) for a discussion of the various frames that can be used 
to describe bubble collisions). This yields the position z = Zc oi the collision. 

2. At the collision, impose energy conservation to determine if a solution for Mq and M^ is possible for a set of 
input parameters: Ha, Hb, He, Zc, Rq, ksc- If a- solution exists, this provides the initial conditions for the 
domain walls enclosing the C phase. 

3. Solve for the motion of the domain walls enclosing the C phase using Einstein's equations. 



Junction condition formalism 



A domain wall is a (2-1-1) dimensional boundary between two vacua, with metric functions a^ to the left and an 
to the right. The geometry is continuous but not smooth, so the z coordinates are the same but xl, xr can be 



different. The normal for the wall is defined to point from the L phase into the R phase, and requiring orthogonality, 
its components are given by 



L 



XL, 



R : 71^ = ifl. 



(4) 



where the dot refers to a derivative with respect to the proper time of an observer riding on the wall. Normalizing, 
we obtain the condition 
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Integrating Einstein's equations across the wall, we obtain [10], 

Klb~K%b = -kLR^t , 
where K is the extrinsic curvature of the wall (which is different on either side) , 

1 



K„ 



'd^-fab , 



(5) 
(6) 
(7) 



kLR = inGaLR. is the normalized tension of the wall, and a, b run through the (2+1) dimensions on the wall worldsheet. 
Although the components of Eq. ([6]) appear to yield two independent equations, they are actually the same and simplify 
to 



Pl- Pr = kLRZ 



(8) 



Here /3/ — aiij is defined as the physical "boost" of the wall, which will simplify further calculations and their 
physical interpretations. 

From Eq. [s] and the definitions Eq. (Is]) , we get 
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2.h. ■ <"> 

In this form, we can already note an interesting property when a^ and ur are positive definite. In this case, the sign 
of XL. R are determined by the sign of /3l. r- When \aR — a^j > z^fc^, we have xl, r > or xl, r. < 0, and the wall 
moves away from the region with smaller a/ into the region with larger a/. We will refer to this behavior as normal. 
When \aR — ulI < z'^k'^, we have xl > while xr < 0, and the wall moves away from both sides. We will refer to 
this as repulsive behavior. 

Integrating Eq. (lo]), it is possible to solve for the trajectory x{z). Note that Eq. ^ is just the equation of motion 
for a particle of unit mass with zero energy moving in an effective potential 
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In general, one can consider initial conditions with either positive or negative i. To better understand the motion in 
the effective potential, we can manipulate it in 2 ways: 

- 8K// = ^ 
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(13) 

(14) 



It can be seen that both a^ < and ur < at some position z are necessary conditions for a turning point in the 
motion of z {z — 0). Otherwise, the effective potential is negative definite, and the motion in z must be monotonic 
from z = Zc — > oo or from z = Zc — >■ 0. 




FIG. 4. Our coordinates only cover the interior of the forward lightcone of the x axis. Therefore, a full S0(3,l) symmetric 
bubble wall shows up as two separated portions, namely two separated lines on the x — z picture. Each point there is a 
suppressed H2 on the left picture. Only the blue(right) portion will be relevant to a collision with another bubble from a larger 
X position. 



C. Determining the kinematics of the collision 



Before the collision, the incoming domain wall in the S0(2,l) coordinates is only a portion of the S0(3,l) symmetric 
bubble wall, as shown in Fig. |4] By symmetry, we can choose the incoming bubble on the left to fully describe the 



initial conditions for the collision. With ol 
the incoming domain wall are quite simple. 



as 
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In order to collide, the incoming wall has to move into region A, which requires xa > 0. The wall trajectory always 
begins at 2: = where a a > 0, and from Eq. Mm we require fc^^ < H\ — H^ so that the wall is normal and Rq < H^^. 
If the bubble centers are initially separated by a distance 2Aa;, then arranging the incoming walls so that they 
collide at X = 0, the solution to Eq. ( 15 1 for the trajectory of the wall is given by 



z = H^^ [(1 - HlRl) sec2 [Ha (a 
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if Ha = 0. The location of the collision, z — Zc, is found by setting x = in these expressions. 



IV. COLLIDING DE SITTER AND MINKOWSKI BUBBLES 



In this section, we look for classical transitions in cases where the bubble interiors have Hg > 0, leaving Hq a 
free parameter. The collision occurs at {x = 0,z = Zc}, and because the colliding bubbles are identical, there is a 
reflection symmetry about x = 0. We can therefore choose to describe the collision by focusing on the incoming and 
outgoing walls on the left side. The motion of the incoming walls was found in Sec. |III C[ After the collision, in the 
absence of null shells of radiation, the post-collision domain wall is described by 
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(21) 
(22) 



In the B region, the metric function Ob is positive definite for Hg > 0. In addition, by symmetry, we must have 
ac > 0. Therefore, in order to have a region of C that is initially growing after the collision {xq < 0), we must 



impose the physical constraint Pc < 0. Since gb is positive definite, we have from Eq. ( 14 1 that after the colhsion z 



must be monotonically increasing (z is a timehke coordinate in B, and therefore initial conditions with i < are not 
physical). 

The parameter Mc in ac = 1 — 2Mc/ z + H'^z^ is determined by energy conservation. Technically, it says the boost 
angles around the collision point Zc, after moving across all crossing walls, must sum up to zero [T51 - fT5] . 
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In the above equation, all functions of z refer to their specific values at the collision point z^. Also, since the incoming 
boosts are replaced by avi/, we will drop the "out" subscript and (3 from now on refers to the outgoing walls. Further 
simplification leads to 



(24) 



Although simulations [SI [9] show that in single field models very little energy is radiated in a classical transition, 
the addition of extra degrees of freedom might allow for significant energy losspU]. A detailed picture of how much 
energy is lost and in what form is beyond the thin-wall analysis that we perform. However, the main impact is to make 
a classical transition harder by taking some energy away. In an attempt to model this, we include shells of radiation 
as shown in FigjSJ This choice of energy sink is motivated by simulations [I H 118 ) . In this case, the additional region 
has the metric function 
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where Mg is related to the energy density of the radiation shells by 
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The boost of outgoing walls takes exactly the forms in Eq. (28) with a^ replaced by a^. 
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The energy conservation equation picks up an additional term. 
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which in the end is very similar to Eq. ( 24 1 
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Note that with the extra variable Af^, there are not enough constraints to fully determine the geometry, and we 
need detailed information about the microphysics as discussed above. In practice, we will treat Mj^ a tunable input 
parameter. 



Null limit 



There is a particularly simple form of the energy conservation equations when the incoming and outgoing domain 
walls are highly boosted. The incoming walls will be highly boosted when they are light {k\g ^ H\ — H^) and z is 
large {zc ^ H^ ). In this case, we have that aw ^ aA, aw ^ aB- The out going walls will be highly boosted when 



—Pc ^ ac and — /3^ ^ a^ (or — /3^ ^ a^ if radiation shells are included). In this limit, both Eq. (24) and Eq. (30) 
simplify to 

aAac = a% . (31) 

This is exactly true for null domain walls and numerically verified to be a good approximation for highly rclativistic 
walls. We can use this simple relation to examine the general properties of classical transitions. 
In this limit, ac{zc) is independent of H^, which is apparent after solving for the value of Mc- 

^ -^^"^'^- l + H\zl ■ ^^^> 

For Hg < H\, we always have ac < as at the collision point, and the domain wall is normal, moving from region 
C into region B. This makes sense as the collision is a local event, which only concerns the incoming and outgoing 
walls, but not the properties of the spacetime regions in between (such as their vacuum energy). 

However, as z increases at late times, the Mc term in ac becomes irrelevant and only the vacuum energy matters. 
In this regime, the signs of xb and xc depend only on the values of the energy density in either vacuum and the 
tension of the post-collision domain wall: 

sign(iB) ^ sign(i/2 -Hl + k%c) , (33) 

sign(xc) ^ sign(H2, - ij| - k^c) ■ (34) 

Thus, adjusting H^ controls the asymptotic behavior of the outgoing walls, and produces a complete list of possible 
geometries as the results of classical transitions from identical bubble collisions: 

• H]j - H^ > kg(j Normal 

After the collision, the domain walls always move from region C into region B {ic < 0, xb < 0) . Namely, they 
are always normal. 

• \Hg — _ff^| < kj^fj Repulsive 

The initially normal domain walls eventually become repulsive. As a result, an expanding region C seems to 
be moving away from region B on both sides of the wall (xc < 0, xb > 0) . 

• H^ — Hg ^ kj^fj Oscillatory 

Immediately after the collision, the domain walls move from region C into region _B, but eventually turn back 
and collide with each other. The second collision can induce another transition from C back to A with similar 
domain wall motion, resulting in an oscillating series of transitions. In Appendix [X] we outline the properties 
of these geometries in more detail, showing that each successive collision is less and less energetic. This loss of 
energy can cause the oscillations to terminate, at which point the collision energy must be dissipated in some 
other form. If the oscillations last long enough, the oscillating region acts as a Repulsive domain wall separating 
the interiors of the two colliding bubbles. 

• H^ — Hg > kgfj Marginally Repulsive 

A while after the collision, the initially expanding domain walls turn back to move from region B into region 
C {xc > 0, Xb > 0)- However, in classical transitions with H^ > 0, the expansion in region C and the 
gravitationally repulsive nature of domain walls prevents the walls from colliding again. 



Finding exact solutions to the energy conservation equation Eq. ( 24 1 , we then solve the equations of motion for x 
and z (Eq. [9]) numerically to confirm the existence of each of these geometries. Such numerical simulations are key to 
determining if and when the Marginally Repulsive and Oscillatory geometries exist, since they depend on the details 
of the evolution away from the collision. Examples with Hj^ > and Hq > are shown in Fig. pj In this series of 
plots, we show the wall trajectory as a function of the conformal time 

HjTj^ I dz— . (35) 

J ai 

In these coordinates, future infinity is located at HjTi — 7r/2. The location of a collision is denoted by the solid disc, 
and the post-collision walls correspond to the red dashed line. The interior of the second colliding bubble is a mirror 
image of the plots on the left. 

The Oscillatory solutions are described in detail in Appendix!^ and here we show the first transition to the region 
of C. The field dynamics suggest that the secondary collision will yield a classical transition back to A. For the type 
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FIG. 5. Examples of the Normal, Repulsive, Oscillatory, and Marginally Repulsive classical transition geometries. The left 
figure is the B bubble on the left side of the collision, with its initial domain wall to vacuum A(blue, solid) and the domain wall 
to vacuum C(red, dotted) after the collision. The right figure is the middle region C with its domain walls(red, dotted) and the 
event horizon(black, solid). In all but the Oscillatory geometry, an infinite region of the C vacuum is produced by the classical 
transition. In each case, we choose the parameters (in units oi Ha) Normal: {_Ro — .01, Aa; = 2.9, He = -SjHb — ■2,k = .01}, 
Repulsive: {Ro = .01, Ax = 2., He = .4, ffs = .2,k = .4}, Oscillatory: {Ro = .01, Aa; = 2., He = .3, Hs = .2, fc = .01}, and 
Marginally Repulsive: {Ro = .01, Ar? = 2.9, He = .3, Hb = .2, fc = .01}. 
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of potentials we consider, the oscillations will eventually terminate, and the bubbles of B will merge. The Marginally 
Repulsive geometries are more interesting for the question of which vacua can be populated in a classical transition. 
In our numerical simulations, the existence of these geometries can be easily verified by first drawing the past light 
cone from the point x = at future infinity in the region of C. If the wall exits this light cone, it can never recollapse 
to zero size. Therefore, even though xc > 0, the region of C will continue to grow in physical volume. We have 
verified the existence of these geometries over a wide range of parameters, and we show a representative example in 
Fig.[5] 

In the absence of gravity, only the Normal and Oscillatory geometries exist, and a lasting region of the C vacuum 
can only be produced if H^ < Hg. Including gravity, the Repulsive and Marginally Repulsive geometries allow for 
lasting regions of C even when H^ > Hgl This is due to a combination of the gravitationally repulsive nature of 
domain walls and the background expansion of the C vacuum, and is a key result of this paper. We now discuss the 
range of parameters in it is possible to find this behavior. 

B. Constraints on possible geometries 

Here we will go beyond the null approximation to investigate parameter constraints on the possible geometries, 
especially for the (marginally) Repulsive geometry. The constraints come from simple physical intuition — energy 
conservation prevents the creation of arbitrarily heavy objects. Although the incoming domain walls can accelerate 
for an unbounded amount of time, the expansion of de Sitter vacuum A still makes the total energy (in the center 
of mass frame) bounded. We should therefore expect a bound on the tension ksc of the post-collision domain walls 
that a collision is able to produce. Let us now see how this arises. 

For simplicity, we still demand the vacuum bubbles of B to have small critical size, Rq <S^ H^^ , therefore the second 



bracket on the right hand side of Eq. (|30|) is approximately 1. Using the relation Pb = Pc + kg^z'^, we obtain 

(36) 



^/aA{a^ + z'^k^.^j + 2zrkBcPc) _ -Pc + \/ac + Pc 



2 



0-B -Pc ~ Zcksc + \/ag + [Pc + Zcksc) 

The physical range that (3c can take goes from f3c = 0, in which case the post-collision domain walls are produced at 
rest, to f3c — — oo, where the walls are highly relativistic. Plotting both sides of this equation as functions of /3c in 
the physically allowed range, f3c < 0, we get Fig|6] The two curves have to cross somewhere in this range in order to 
have a solution. As we can see, the left hand side goes to zero at some finite (3c while the right hand side asymptotes 
to 1 a.s (3c —>■ — oo. Therefore, the existence of a solution implies that at (3c = 0, the left hand side is larger than the 
right hand side. 



as -Zcksc + V°^g~+^P|c 
After several lines of algebra, we get 



"' > . . ' , „., ^ m 



kBc < J" r^ ■ (38) 



This agrees perfectly with our intuition: 

• Setting (3c to zero saturates an upper bound on ksc- That is, the heaviest possible out going walls are produced 
at rest. 

• Energy leaks into shells of radiation imply a positive Mg, therefore a smaller a^ < as, which makes the upper 
bound for ksc smaller. The addition of energy sinks makes classical transitions harder. 

• As Zc increases, the bound on the tension gets larger but asymptotes to 

Max(fcBc) = ^^1^ ■ (39) 

The center of mass energy is bounded for H\ > 0. 

Now that we have specific constraints on the possible classical transition geometries, we can return to the question 
of when a classical transition can produce a lasting region of vacuum C. Again, we are particularly interested in cases 
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FIG. 6. The blue curve is left hand side of Eq. (36 1, which goes to zero at some finite /3c- The red curve is the right hand side, 



which asymptotes to 1 as Pc -^ — oo. The parameters here are Ha ~ 2, Hb = 1, He = 0.2, fc^c ~ 0.4, z^ = 10, Ro = 0.1, 



Mg = 0. 



where He > Hb, since gravitational effects will be important in this regime. One way to guarantee a surviving region 
of C is to produce a Repulsive geometry, which requires that ksc > Hq — Hg. Combining this with the maximum 



tension Eq. (39 1, we obtain a bound on the vacuum energy of the C: 

v2 



Hh < ^^^^^ (40) 

For vacuum energies satisfying this bound, energy conservation will allow a Repulsive geometry. Note that by default 
H^ < H\ (which is required for bubbles expand and collide), so this also implies H^ < H\. Numerical study shows 
that the Marginally Repulsive geometry also can only exist for H^ < H\. Therefore, it is possible to produce a 
classical transition to lasting regions of C with H^ > H^ only when H^ < H\. 

C. Embedding regions of higher energy density 

Naively, it might seem surprising that it is possible to find geometries where a lasting region of C can be produced 
with H^ > H^. Viewed from the perspective of vacuum B, this is a false vacuum (high energy density) "blob." In 
spherically symmetric geometries, there is a theorem due to Vachaspati and Trodden (VT) [3T] showing that a false 
vacuum region must be larger than the true vacuum horizon size. At first, our solutions might seem to be in conflict 
with this theorem: the false vacuum region C grows from zero size, and then inflates indefinitely. It is not an exact 
logical conflict as our geometry lacks spherical symmetry, but that is a meek point. Spherical symmetry was only a 
technical convenience used to bypass the complication of caustics in the theorem. We would like to provide a better 
physical reason to reconcile the apparent paradox. 

Let us examine the assumptions of the VT theorem. First, the Null Energy Condition (NEC) must be satisfied. 
VT additionally assume spherical symmetry, and then the Raychaudhuri equation for the divergence of a spherical 
congruence of null geodesies can be written as 

dO , , 

The content of this equation is that initially converging in-going null rays cannot diverge unless they pass through 
the origin. Said differently, following a null ray, the surfaces it passes through cannot go from being normal to 
anti-trapped. For a spherically-symmetric region of false vacuum embedded in some background with H^ > 0, this 
condition implies that the initial size of the false vacuum region must be larger than H^^ . 

The bubble collision geometry we are considering is shown in Fig [7J which depicts a constant time slice through 
the classical transition. Consider the case where the tensions are small, and the vacuum energies are comparable 
in magnitude. We can then define a set of flat coordinates, approximately given by the flat slicing of the exterior 
de Sitter space, whose origin is located at the center of the region of C. In these coordinates, the physical volume 
removed from the B phase is equal to the physical volume replaced by the C phase. Now, draw a closed spherical 
surface, whose radius is smaller than H^^, but larger than H^^, that completely encloses the region of C at some 
time. Following a congruence of null rays inward from this surface, rays in the B phase will be converging, while rays 
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Hc<R<Hb' 



FIG. 7. The classical transition geometry on a constant-time slice. The region in the C phase is shaded (blue). A surface of 
1 ^ D ^ Tj- ^"^^i encloses the region of C is shown as the solid circle. A spherical congruence of ingoing null rays 



radius H(^ < R < Hg 



an instant before and after this time slice are depicted as the dot-dashed and dashed lines respectively. Diverging portions of 
the congruence are denoted by the solid red lines. 



in the C phase will be diverging. The shear and twist induced across the various interfaces wiU be minimal as long 
as the vacuum energies are comparable and tensions small. Now, tracing a congruence of null rays that move inward 
towards this surface, null rays in the B phase will again be converging, but as long as Ha > He, the null rays in 
the A phase will be diverging. The surfaces that the congruence of null rays pass through can therefore always be 
anti-trapped as long as Ha > He, and there is no conflict with the VT theorem. The reason it is possible to embed 
a region that is higher energy than the B phase is that in the collision geometries, the A and C phases are always in 
contact. 

There is another, related, theorem put forward by Guth and Farhi |22| : embedding a small and growing region of 
false vacuum in asymptotic Minkowski space must be accompanied by an initial singularity. This is an application 
of a singularity theorem due to Penrose [23; , and does not require spherical symmetry. It does however require the 
existence of a noncompact Cauchy surface — a condition our geometry naively violates by having a de Sitter parent 
vacuum with Hj^ > 0. However, Farhi and Guth made a reasonable conjecture, based on the fact that local physics 
should not care about the cosmological scale HJ^^. They claimed that the same theorem holds for de Sitter spaces 
as well, if all our operations to generate the false vacuum region are "local" . Indeed, it can be demonstrated that 
regions of false vacuum embedded in a background de Sitter space either require the existence of an initial singularity, 
or must be larger than the exterior horizon size |24j . Our geometry pushes the conjecture of Guth and Farhi to the 
extreme, and verifies their claim. Since the colliding bubbles can be almost 27?^ apart, the initial configuration of 
the two colliding bubbles can be considered "local" only in the sense that the bubble walls eventually collide, and 
produce a classical transition, in the asymptotic future. 



CRUNCHING BUBBLES 



In this section, we focus on the case where H\ > and H^ < 0. The interior of bubbles with H^ < undergo 
a big-crunch [25j . which has earned minima with negative energy the colorful title of "terminal" vacua - this is 
supposedly the end of the line for any observers inside of the bubble. But is it? We have seen that classical transitions 
can produce a lasting, non-singular, inflating region of C. Could a classical transition allow some portion of the 
bubble interior to avoid the crunch? At first glance, the prognosis is encouraging. For example, by taking Ha = 0, we 



can have arbitrarily large center of mass energy in the collision, and the bound on the tension ksc in Eq. ( 39 1 does 
not exist. This implies that it might always be possible to create heavy walls, and therefore Repulsive geometries, in 
which case the region of C, regardless of its vacuum energy, is guaranteed to survive. 

However, producing an inflating region of C from the collision of bubbles in Minkowski space is in conflict with 
< H\ from the previous section. One technical twist in the case of AdS bubbles is that there are 



our bound H^ 
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FIG. 8. The Penrose diagrams for AdS(left) and Hyperbolic-Schwarschild (right), for Mc > 0. The x and z coordinates are 
increasing in the directions labeled by the arrows. The blue curve (left-moving in the AdS diagram) has i > and the red 
curve (right-moving) has i < 0, but they have the same /3b ~ asXB > 0. 



horizons, across which z becomes the spacelike coordinate and x becomes timelike. See FiglSlfor a complete coordinate 

cover. For z < Hl^ 

-1 



'B 



the matching problem is exactly the same as Hp > 0. 



LB — ^- Highly boosted domain walls will collide 
in the region where Zc > H^'^. Something interesting happens there, and we will see that the bound Hq < H\ 
on surviving regions of C is enforced not through kinematics, but through the large-scale structure of the classical 
transition spacetime. 

In Sec. |V A we will first demonstrate how the different matchings work and how they depend on various parameters. 
Then in Sec. VB we will discuss the consequences: It is possible to make a lasting inflating region C in the collision 
of two AdS bubbles, allowing some regions to escape from the crunch. However, such regions cannot have a higher 
vacuum energy than the parent vacuum A, just as we found before in the collision of de Sitter and Minkowski bubbles. 



A. General Matchings 



We start from an example with H^ = 0, which will capture most of the physics we want to understand. The 
formulas we derive, however, will be completely general. The metric in region C is Hyperbolic-Schwarschild, with 
metric function ac = 1 — ^^7^ • The causal structure is drawn in Fig. M 

In the region where z^H"^ > 1, (3b is guaranteed to be positive, since as < and a; is a monotonically decreasing 
timelike coordinate. There will be two timelike vectors with the same Pb , one with positive z the other negative. This 
is not an ambiguity, but really two physically distinct situations. When i > 0, the collision in region C is located in 
the expanding quadrant; when 2: < 0, the collision is in the shrinking quadrant. 

In terms of equations, the i > case is simpler and is where we shall start. The analysis is nearly identical to 
the dS and Minkowski matchings we studied earlier, but for convenience we will rearrange the energy conservation 
Eq. 1241) to 



y/OAOC 



-as 



law + ^aw — o-A 



aw + \Jaw -as) \ I3b ~ \//3| + ub 



-Pc + Vl^c + ac 



(42) 



so that both sides are positive definite. Unlike the collision between dS and Minkowski bubbles, there is no null limit 
of this equation where /3c — >■ — cxd. This is due to the fact that xb is a timelike coordinate for Zc > \Hb\~^, which 
from Fig. [8j implies that /3b > 0. The condition /3b — (3c + fcsc^ then shows that (3c cannot become arbitrarily 
negative. 

Following the analysis of Eq. (36), by looking at Eq. (42 1 in the physical range /Jp < and ac = 0, we can derive 



bounds on kBc- Although the starting equation looks identical, the new roll of xb as a timelike coordinate yields an 
interesting twist: comparing the two sides of the equation, in addition to the upper bound, we also get a lower bound 
on fcsc: 



UA + (-Qb) 



> ZckBC > 



-as 

y/oA 



(43) 
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For this range to exist, we need (—as) < cla- Note that in order to guarantee xb timehke, ksc is aheady bounded 

^I^'bc > ac-a-B > -aB ■ (44) 



Under the condition that (—as) < ga, this lower bound is strictly better than Eq. (43), so logically we should replace 
it. 

> zJiBC > Va-c - as ■ (45) 



2y/aA 

This lower bound depends on a dynamical variable ac, therefore cannot truly be part of a physical constraint. It will 
soon become obvious that this bound comes from our assumption of i > 0, while nothing spectacular happens when 
z changes sign. 

For i < 0, as can be seen in Fig|8J the solution should have /3c > because xc is increasing as one moves out of 
the interior of the region of C. The corresponding equation energy conservation is in this case slightly different from 



Eq. (24) (the first factor from incoming walls is already set to one). 

yaAac _ \/Pc + ac + ^c 
-as V^l + as + Pb 



(46) 



The first good sign here is that the null limit exists, and it is exactly aAO-c = a^ as before. For a physical constraint 
on kBCi we again compare both sides of Eq. ([46]) as functions of j3c over its physical range Pc > 0. When /3c — > oo. 



the right hand side of Eq. (46) goes to one, while the left hand side goes to infinity. A solution is guaranteed if 



there exists some value of /3c > where the right hand side is greater than the left hand side. The strictest bound 



on the tension comes when Eq. (46) is evaluated at some physical lower-bound on f3c- There is always a solution 
if kgQz'^ < (—as), since the physical lower bound is oc = in this case. No lower bound on the tension exists for 
solutions with i < 0. If ac — as > kg^z'^ > —ub, then the physical lower bound on /3c is /3c — 0, which leads us to 
the same upper bound on /c^c s-s we obtained for the i > case, namely: 



aA+ -as) 
ZckBC < ?r^ — , (47) 



knc < ^±So^ , (48) 

ZZcJaA 



but with the requirement a a < {—clb)- 
Summarizing the results of this analysis: 

• The kinematic upper bound is always 



which takes exactly the same form as in Eq. ( |38[ ) r] In Minkowski space, this upper bound does indeed become 
arbitrarily large. 

• Although matchings with different signs of i looks quite different in FigjSJ they are continuously related situa- 
tions. For light domain walls, kBC'Zc < \/—o-b, we always get i < 0. For heavier ones, kBc^c > ^J—clb, we get 
i < when a a < (— os), and i > when a a > (— as)- 

B. Escaping the crunch 

Now that we understand how the matching works, we turn our attention to two interesting questions. 

• Can we escape from the AdS crunch in the colliding bubbles by creating a classical transition to a lasting region 
with H^>07 

• Can the AdS crunch play a role, much like the initial singularity was important for 'creating a universe in the 
lab' [22], possibly allowing for solutions with Hq > H\l 



Radiation leaks have been ignored in the AdS analysis but it is easy to see that restoring them reproduces Eq. MSk 



16 

Recall that the way to guarantee a lasting region of some vacuum C is to create heavy domain walls, and hence 
Repulsive geometries. When a^ < {—as), the heavy walls will have i < 0. In this case, the region of C will either end 
in a cosmological crunch, or result in an Oscillatory geometry regardless of the value of the tension kgc- We prove 
this claim in Appendix O In this way, the large-scale structure of spacetime steps in to enforce limits on the types of 
classical transition geometries that can be produced - Repulsive geometries simply don't exist, no matter what the 
tension is. 

The story is potentially different when a a > (—as), where the heavy walls will have z > 0. In this case the 
necessary condition to form a lasting region of C becomes the same as in the de Sitter or Minkowski collisions — we 
need a Repulsive or Marginally Repulsive geometry. Again we prove that the Repulsive case is kinematically allowed 
analytically, and leave the Marginally Repulsive case to numerics. On e key ingredient here is a large tension, fc^c, 
which indeed can become unbounded if H^ = as shown in Eq. (48). On the other hand, the AdS region B also 
increases the threshold for the Repulsive geometry, 

klc>Hh + \HB\^ ■ (49) 



Combining this with Eq. (48), we get 



^ < 4(1 + zlHl) • ('") 

It is easy to see that H^ > is quite achievable, actually quite natural for large Zc- So our answer to the first question 
is yes! A classical transition provides a way to escape from the AdS crunch, at least for some observers. 

For the second question, we further demand that region C has higher energy density than the false vacuum A phase: 

^2 .„2 ^ zl{Hl-\HB\'f~mB? ,.^. 



A few lines of algebra leads us to 



(-as) > ioA . (52) 



Note that {—qb) > aA is enough to guarantee a matching with i < 0, which cannot produce a lasting region. Eq. (52) 
exceeds that threshold by a factor of 3, which strongly suggests that a Marginally Repulsive geometry will not work 
either, a statement supported by numerics. Therefore, the structure of the spacetime never allows a lasting region of 
C to be produced when 7?^ > H\. 

VI. DISCUSSION AND CONCLUSIONS 

Classical transitions arising from the collision between vacuum bubbles provide an alternative way of populating 
regions of the landscape. In this paper, we have found that gravitational effects allow for classical transitions to 
populate a much wider variety of vacua than simple field dynamics suggests. In particular, it is possible to produce 
regions of vacuum with a higher energy density than that of the colliding bubbles. Even collisions between crunching 
bubbles can produce a lasting region of inflating vacuum via a classical transition. What could these findings mean 
for populating vacua in eternal inflation? 

One important aspect of classical transitions is that vacua which cannot be accessed by bubble nucleation can 
be reached through a classical transition. For example, considering only bubble nucleation, an AdS minimum that 
separates two sets of inflating vacua (in a one-dimensional potential) can carve the potential landscape into two 
disconnected islands [2S] . Classical transitions can provide a bridge between such islands, allowing a wider variety of 
vacua to be populated during eternal inflation. 

Even for vacua that can be reached through bubble nucleation, classical transitions might be the dominant pro- 
duction mechanism. This could occur if the bubble nucleation rate directly to the vacuum of interest is quite slow, 
but the rate to an adjacent minimum is relatively fast. Collisions between the quickly-nucleating bubbles might pro- 
duce transitions to the slowly-nucleating vacua much more efficiently than direct bubble nucleation. Estimating the 
probability for a classical transition to occur as the joint probability of two bubbles nucleating in the same Hubble 
patch within the same Hubble time (so that they collide) , this occurs roughly when the nucleation rates per unit four 
volume satisfy T\g > Vac- 

In both of these cases, classical transitions have the potential to play an important role in the measure question for 
eternal inflation. The starting point for most measure proposals is a matrix of transition rates between the vacua in a 
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landscape, which is then used to determine a regulated distribution of volume, bubble number, or some other quantity 
that can ultimately be tied to the relative number of observers in different vacua (see e.g. [H [5] for a review). Altered 
transition rates notwithstanding, depending on the choice of measure, classical transitions can drastically affect the 
overall picture of eternal inflation. For example, if one chooses a particular bubble nucleated during eternal inflation, 
all but a set of measure zero of the volume on a constant FRW time slice inside the bubble will be to the future of a 
collision [27j . If collisions can produce a classical transition to a lasting region of some new vacuum, then effectively 
all of the volume on any constant time inside the bubble will be removed, and replaced by the new vacuum. This is 
true both in the case of terminal (AdS or stable Minkowski) and non-terminal (de Sitter and unstable Minkowski) 
vacua. We hope to explore the implications of these and other issues in future work. 

Regarding the "no universe in a lab" conjecture [2TJ [22], our geometries provide some new insights. It is not easy 
to find explicit geometries without spherical symmetry, which are easily calculable, and support the conjecture. The 
fact that classical transitions with He > Ha cannot occur from the collision of two bubbles, a fairly non-local and 
non-spherically symmetric initial condition, should dramatically boost our confidence on the conjecture. 

Another area for future work regards the Oscillatory geometry. In our analysis, we have found that energy con- 
servation allows for many oscillations, and at late-times the region undergoing the classical transition acts as an 
effectively Repulsive domain walQ This could lead to interesting observational signatures of the collision, since there 
is a somewhat stable topological defect that can persist during inflation. 
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Appendix A: Details of the Oscillatory solutions 

In this appendix, we treat the properties of the Oscillatory solutions in more detail. As we discussed in Sec. |IV[ the 
initial collision between bubbles of the B phase can give rise to a classical transition to the C phase. In the Oscillatory 
geometries, the post-collision domain walls enclosing the C phase eventually re-collide. At the second collision, the 
fleld dynamics suggests that if a classical transition occurs, it will be to the A phase. This region can expand and 
re-collapse, possibly giving rise to a classical transition back to the C phase, and so on, producing a pattern of 
classical transitions C — ^ A — >■ C — >■ ^ . . . . Imposing energy conservation at each collision, we can numerically evolve 
the equations of motion on each segment, and determine the full solution. The main conclusions we draw from this 
set of numerical simulations are 

• If the walls separating B from C and B from A have different tensions, there will be a finite number of oscillations. 
Eventually, there will not be enough center of mass energy in a collision to produce the heavier set of walls. 

• When the oscillations terminate, there are two possibilities. First, there could be continued oscillation between 
the same vacuum: C — A — C — A — A— ...or C — A~C~C — ... depending on which set of walls is heavier. 
Second, there could be energy loss through an ejected shell of radiation. In this case, eventually the two colliding 
bubbles merge. 

• If the walls can oscillate for long enough, the region of C and the walls that enclose it become one effectively 
repulsive domain wall: that is, xb < 0, and the walls are always moving out of the B phase in both of the 
colliding bubbles. 

To begin, we index the collisions by z = 1, 2, 3, . . ., where i = 1 is the initial collision between the bubbles of B, and 



j > 2 label the subsequent collisions of the Oscillatory geometry. We must impose energy conservation Eq. ( 24 ) at 



This has been explored in the braneworld scenerio with 1 higher dimension|28| 
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FIG. 9. The outgoing velocity xc (left) and xb (right) for a number of collisions i = 1, 2, ... 19 in the Oscillatory geometry. This 
simulation was for the parameters (in units where Ha ~ 1) {^o = .33, Ax = 2.2, He = A,Hb = .01, /cas = .05, /sab = .075}. 
Note that the envelope of xc decreases with i. The classical transitions must either terminate, or only involve vacuum C for 
i > 19. For j > 7, xb > 0, and the walls enclosing the region of C move out the interior of both bubbles of B. 



each collision, which for i > 2, and neglecting the presence of any other energy sinks, is given by 



^ ac-iflc. /-fc-i + y^/?c.-i + «c.-i \ /-/3c. + y^/^c. + 



as 



-/3b._, + ^P%_, + 



a-Bi 



«c. 



-/3b. + V/3|^ + OB. 



(Al) 



Each of the quantities are evaluated at the position of the i-th collision z — z^ (by symmetry, all collisions will occur 
at Xc = 0). The metric coefficient a^ generically denotes the region undergoing the classical transition due to the 
i-th collision, and is given by 



aci 



1 



2M; 



Hlz\ 



(for i odd), ac;^ = 1 



2M. 

z 



H\z\ 



(for i even) 



(A2) 



To satisfy energy conservation, the mass parameter Mi must be different for each classical transition. Conservation 
of energy at the «-th collision then provides an equation for Mi in terms of Mi_i. The definitions of /3b. and Pd are 
as in Eq. (22 ), with ac ^^ a-d ■ 



To find the Oscillatory geometries, we solve for energy conservation, numerically evolve the equations of motion for 
the wall to find the location of the following collision, and iterate. The first observation is that the outgoing velocity of 
the wall, as measured by xc, decreases each time a collision occurs (hence /3c also continually decreases). This implies 
that the collisions become more frequent with increasing z, and that each subsequent classical transition produces 
a region whose maximum size along the xc direction is continually decreasing. Recalling that /3b = /3c + k^z^, as 
/3c (which must be negative) decreases in magnitude, eventually /3b must become positive. When this occurs, the 
walls move away from the interior of both the B bubbles, just as for the Repulsive geometries described in Sec. |IV A[ 
The region undergoing the classical transition acts as an effective domain wall between the two bubbles, which at 
late times has a Repulsive behavior. In addition, because the outgoing velocity is continually decreasing, the effective 
center of mass energy available in each collision is decreasing. Therefore, if the tensions kAB and kBC are different, 
then the total number of oscillations must be finite. Eventually, the center of mass energy will become so small that 
it will be impossible to produce the heavier set of walls. 

All of these behaviors are observed in numerical simulations. An example is shown in Fig. [9] where we plot 
the outgoing value of xb and xc for a number of collisions. Here, the tensions kAB and /cbc are not equal, with 
kBC < kAB- The red points correspond to classical transitions to the C phase (for the collision index i odd), while 
the green dots are for classical transitions to the A phase (for the collision index i even). It can be seen that xc 
oscillates, due to the fact that fcsc < ^ab, but the envelope continuously decreases in magnitude. For the collision 
i = 18, Xc is nearly zero, and there is no solution to the equations for energy conservation at i = 20. The classical 
transitions must either terminate at this point, or transitions must occur to the C phase only. Studying the plot for 
Xb, it can be seen that for i > 7, xb is positive. Therefore, for the duration of the subsequent oscillations, the walls 
are continuously moving away from the interior of both bubbles of B. 

Interestingly, in some potentials the Oscillatory solution must last forever. The simplest example is a special case 
of the model presented in (llj , a single scalar field with a Z2 symmetric potential that has a false vacuum A in the 
middle, Bi to the left and B2 to the right. The Z2 symmetry implies that Ha > Hbi — Hb2 and kABi = kAB2- The 
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classical transition in this special case will produce an oscillating region of the A vacuum. Because Bi and B2 are 
different vacua, there must always be something separating the interiors of the two colliding bubbles. The classical 
transition can play this role if the oscillation does not end, and gives rise to a stable configuration in field space that 
goes between the position of vacuum Bl, through A, to the position of vacuum B2. As above, such a configuration 
will move out of both colliding bubbles at late times, producing a Repulsive geometry. 

Appendix B: Collisions between different bubbles 

In this appendix, we consider the case where the colliding bubbles that produce a classical transition are different. 



The equation of motion for the pre-coUision and post-collision domain walls are identical to those presented in Sec. IV 
However, the equations of motion will be different for each of the four walls, and we must evolve them separately. The 
condition for energy conservation will also be different. We suppress the full expression, but it can be obtained as a 
straightforward generalization of the analysis in Refs. P"5HT5] . In the limit of null walls, energy conservation becomes 

aAac = aBiaB2 , (Bl) 

where 

aBi = l + Hl,z^, aB2 = l + Hl2z\ (B2) 

are the metric functions in each of the two colliding bubbles (where Bl denotes the bubble nucleated at x < and 
B2 the bubble nucleated at x > 0), with a a and ac defined as before. Also as before, we can define the functions 

ac - aBi - ^^fc|ic « aB2-ac + z^k%^c (T,-i\ 

Pci = 7)—, ' PC2 = TT-, , [D-i) 

2zkBic 2zkB2C 

where the physical criteria for a classical transition to occur is 

fic2 - /?ci > . (B4) 



Substituting for ac using Eq. (Bll, we have in the limit of large Zc that 



^ ' '^2 , . tt2 -"B2 \ ^ ^ Tj2 ^B2 tj2 ,2 



For H\ > H%^ g2i which we require to have vacuum bubble solutions, the left hand side is positive definite, while 
the right hand side is negative definite. Therefore, in the null limit there can always be a classical transition. 

The classes of geometries that can be produced are similar to those produced in the collision between identical 
bubbles. 



• 



If both walls are accelerating away from C, (which includes being Repulsive or accelerating into Bi), then region 
C is lasting. 



• If Hq — Hg^ ^ ^Bic then the geometry will be Oscillatory. However, in this case, the subsequent collisions 
will not occur at a; = 0, but will initially drift towards the vacuum of higher energy. This can be understood by 
looking at the z-position of the turning point in x for both walls: 

(B6) 



• 



''1 tt2 _ tt2 _ l2 ' -"2 tt2 _ tt2 _ U2 

"C -"bi '^BIC ^C "B2 '^B2C 

To the extent that the tensions can be ignored, it can be seen that the turn-around on the side with a lower 
H^-^ will happen first. This means that the wall initially moves towards the bubble containing higher energy 
density. Eventually, just as for identical bubbles, we will have xb > 0, and the wall will accelerate out of both 
bubbles. Also, just as for identical bubbles, since the tensions on either side of the region that has undergone 
the classical transition are not equal, the oscillations must terminate. 

When the two walls behave differently, the general fate of region C depends on the integral 



PC2 - Pc 



ac 



- dz . (B7) 



If it ever goes to zero, the two walls recollidc. If not, it is a lasting region. Nevertheless there will still be a 
general drift into the region with higher -ff^j. 
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FIG. 10. The two triangles in the middle are the full coordinate coverage for Hyperbolic-Schwarzschild with M < 0. The i < 
matching to AdS bubbles B on both sides has to be in the lower triangle, which ends up in a crunch in the same time as the 
AdS crunch outside. 



Appendix C: General geometries with i < matchings 



In this appendix we will show that i < matchings can never have a lasting region C. As shown in Sec jIV A{ the 
null limit of the energy conservation equation is a useful tool to find qualitatively different geometries. 



2Mr 



H. 



2 „2 



\Hb?z^^^ 



1 



HaZ. 



2 ^2 



(CI) 



The parameter Mc can be either positive or negative, and be completely consistent with energy conservation. In 
fact, for i < 0, since we require j5c > 0, the mass parameter must be negative for heavy domain walls. Geometries 
with Mc > are similar to the Hyperbolic-Schwarzschild in Fig. [8J while the Mq < ones are usually drawn as an 
up-side-down triangle with a future singularity. In fact it is equally legal to have a triangle with a past singularity, 
as in Fig. |10[ and a matching with z < corresponds to a portion in this region. In these kinds of matchings, one 
always ends up with a cosmological crunch in the region of C, similar to what occurs inside of AdS bubbles. 

As for Mc > 0, the geometries fall into two classes. Let Zc be the horiz on in the C region (note, this is not 
the cosmological event horizon, but the z-position where ac = 0), Eq. (CI) shows that both Zc > \Hb\^^ and 
Zc < \Hb\~^ are possible. They correspond to slightly different geometries. 

When Zc < \Hb\^^, the domain wall crosses the horizon in region B first. After that ac is positive and decreasing, 
Ub is positive and increasing, eventually we reach a 2 > Zc where 



ac{z) = 03(2) + z k 



BC 



(C2) 



By Eq. (22), this means (3c has to change sign here. The domain walls enclosing the region of C turn around before 



crossing the horizon in region C, and are doomed to re-collide, forming an Oscillatory geometry. 

When Zc > \Hb\^^, the domain wall first crosses the C horizon to where both ac < and ub < 0. Now instead 
of f3c, Eq. (12) shows that i turns around at 2 > \Hb\^^ for which 



aB{z) = {\/-aciz) + zkBcY 



(C3) 



Clearly this is guaranteed to happen before crossing the B horizon. Eq. ( 12 ) also implies that the trajectory is time 



reversible. So after i turns around, it will later enter the top quadrant of the Hyperbolic-Schwarzschild geometry and 
collide with the other wall just like how they emerged initially. This is also an Oscillatory geometry. 
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